Giant Octupole Resonance Simulation 
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Using a pseudo-particle technique we simulate large-amplitude isoscalar giant octupole excitations 
in a finite nuclear system. Dependent on the initial conditions we observe either clear octupole 
modes or over-damped octupole modes which decay immediately into quadrupole ones. This shows 
clearly a behavior beyond linear response. We propose that octupole modes might be observed in 
central collisions of heavy ions. 
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Giant resonances regain much attention presently for 
the investigation of many particle effects in finite quan- 
tum systems. While most of the theoretical treatments 
of oscillations rely on the linear response method or RPA 
methods, large amplitude oscillations require methods 
beyond. Especially the question of the appearance of 
chaos is recently investigated The hypothesis was 

established that the octupole mode is overdamped due 
to negative curved surface and consequently additional 
chaotic damping . Here we want to discuss at which 
conditions one might observe octupole modes at least in 
Vlasov - simulation of giant resonances. 

We will consider different initial conditions of isoscalar 
giant resonances using a pseudo-particle simulation of 
Vlasov kinetic equation . With a local potential U (r) 
the quasi-classical Vlasov equation reads 

-/(p,r,f) + 5_V r /--V p /V r [/(r) = 0. (l) 

at m n 

We represent the distribution function /(p, r, t) by a sum 
of pseudo-particle distributions 

AN 

/(p,r,t)«/ (p,r,t) = ^2-fs(p-p i (t),r-v i (t)) (2) 

i=l 

and use Gaussian pseudo-particles 

f 3 (p - Pi, r - n) = C e - (p - pi)2/2 ^ e -(r-ri) 2 /2-, 2 (3) 

at ri with momentum pi S. These pseudo-particles 
follow classical Hamilton equations 



frp t = -Vf7, r l 



(4) 



We assume for the interacting nucleons a phenomeno- 
logical density dependent Skyrme interaction || which 
results into the mean field 



«7(^) = a— + M — 
Qo \Qo 



(5) 



with a = -356MeV, b = 303MeV and s = 7/6. The com- 
pression modulus is K = 200 MeV. The evolution given 
by (Q) is deterministic, fluctuations appear only due to 
numerical noise M]. We are using 75 pseudo-particles 



per nucleon and a pseudo-particle width, ay = 0.53 fm, 
is adjusted so that the isovector giant dipole energy is 
reproduced at a single mass number. The experimental 
behavior of centroid energy with mass number is than 
reproduced. We have checked different numbers of test 
particles. The dependence of observables on the width 
is discussed in |ll|. Numerically the ground state of nu- 
cleons is realized by Wood- Saxon shapes of density and 
Fermi spheres in momentum. 

The distribution n a (p) = j draif(p,r) of mass distri- 

2 

bution, dj = 1, isospin, etj = r^, kinetic energy, cij = J^-, 
and kinetic isospin energy, ai = -^f, is 
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(2nY 6(p- Pi ) 
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PI 



5(cp - tpi) 



6(0 ~ fr) 
sin(^i) 
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Radial integration determines a spherical distribution 
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f N sinii/i) 
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which can be decomposed into spherical harmonics 

oo / 

Hail), ip) = E E ai m Yi m {d,ip) , 
1=0 m=— I 
AN ^ 
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The observable distributions n a (t9, ip) are normalized to 
y/An aoo, i.e. to mass number A, total isospin T, kinetic 
energy Eki n and kinetic isospin energy EkinT, respec- 
tively. 

As a measure for the strength of the resonances we use 



n„(0)= / dpn«(0,p) = VW^^tcosi?) (10) 

Jo 1=0 
with Legendre polynomials Pi . The amplitudes of multi- 



pole moment, aio y 2l ^ 1 , will be displayed as a function 
of time in the following graphs. The value a 10 means the 
dipole moment, vanishing for isoscalar resonances, 0,20 
is characterizing the quadrupole oscillations and 0,30 the 
octupole ones. 
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As a first initial condition we use the ground state dis- 
tribution of coordinates while the momentum distribu- 
tion is deformed anisotropically. We have modified the 
momenta in a way which corresponds to a giant octupole 
mode. The local densities and currents remain the same 
as in ground state. Figure [ljshows that at start time t = 
there is a pure giant octupole, which is damped out and a 
quadrupole resonance developes instead. The monopole 
and dipole amplitudes which should remain constant doc- 
ument the stability of simulation. 
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FIG. 1. The p icture shows the time evolution of multipole 
moments \f 21 ^ aio(t),l = 0, 1, 2, 3 of mass density concerning 
to normalization, dipole, quadrupole and octupole oscillation. 
Excitation is due to modification in momentum space only. 

In agreement with the already mentioned hypothesis, 
the octupole mode is over-damped. The figure shows the 
nonlinear behavior of mode coupling. Within the linear 
response the damping rate is expected to be independent 
of initial conditions. We choose now other initial con- 
ditions to show that the result is very much dependent 
on initial conditions. Therefore we split the nucleus in 
two parts of mass ratio 3:7 in accordance with symme- 
try of the octupole oscillation and accelerate both pieces 
towards each other. Experimentally it might be realized 
as a central collision of two nuclei with corresponding 
masses. In figure ^| a clear quadrupole resonance appears 
and also a smaller octupole resonance can be seen. Both 
are damped out. Consequently there is no evidence for 
an over-damped octupole mode in this case. 
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FIG. 2. This p icture shows the time evolution of multipole 
moments \/^£r<ilo(t),l = 0,1,2,3 of mass density. Excita- 
tion is due to asymmetric splitting (0.3 to 0.7) of the nucleus 
in two pieces using a plane in coordinate space. Then both 
pieces are accelerated towards each other. 

In order to understand the different initializations, we 
split the kinetic energy into a thermal part and a collec- 
tive part according to (p 2 ) = ((p - (p)) 2 ) + (p) 2 . We 
analyze the time development of the total collective en- 
ergy in the system 



h 2 f 

Ecoii(t) = drg(r)(p) 2 (r,t) 



with the mean current 
1 



<p)M) = 

and the density 
g(r,t) 



dp 



P f(p,r,t) 



dp 



/(p,r,t) 
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(12) 
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In figure || the development of collective energy can be 
seen. There is a background of about 50 MeV due to fixed 
correlations caused by finite width of pseudo-particles as 
one can see from the following estimation. Using (|^) and 



in (11) one obtains 
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x / drf s (r - n(t),a r )f s (r ~ rj(t),a r ) 
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To E E -mPiPifs(n(t) - r,-(f), V2a r ) 
i=i j=i 
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For simulation parameter of 208 Pb, po — 0.162 fm~ , 
N = 75, ay = 0.53 fm, we obtain a basic collective energy 
of 54MeV. Using more testparticles would diminish this 
level. 

The solid line corresponding to figure |l| shows no ini- 
tial collective energy. The exclusive initial excitation in 
momentum space without correlation in spatial domain 
leads to zero initial collective energy. This correlations 
are forming during time evolution. Of course, there is no 
center of mass motion, otherwise we would see just the 
mean streaming velocity. 

This situation is changed if we use the second prepa- 
ration with simple momentum-space-correlations. The 
long dashed line in figure ^ shows initial collective corre- 
lations corresponding to figure ||. Since we can deposit 
enough collective energy in this case, we observe a clear 
octupolc motion. 



cillations where linear response fails. The corresponding 
radius elongation in coordinate space varies about 10 %. 

To summarize, we have observed octupole resonances 
in finite nuclei dependent on the initial configuration. 
The appearance of an octupole mode was shown to be 
possible by correlating the spatial and momentum ini- 
tial excitation. It is possible to excite an octupole mode 
with sufficient collective energy deposited initially. We 
suggest that isoscalar giant octupole resonances should 
be possible to observe in nuclear collisions of mass ratio 
about 3:7 corresponding to octupole symmetry. For a 
mass ratio of e.g. 1:1 we observe no octupole resonance 
in simulation. 

The fruitful discussions with Peter Schuck and P. Li- 
pavsky are gratefully acknowledged. This work was sup- 
ported by a grant (R.W.) and project No. 44185 from 
the Max-Planck-Gesellschaft. 
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FIG. 3. The time evolution of collective energy. The solid 
line corresponds to the excitation scheme of figure |l| and the 
long dashed line corresponds to the figure ti, respectively. 
While the first starts without collective energy, the second 
one starts with maximal collective energy. 

In order to compare with § we calculate the adia- 
baticity index rj defined in j5| as the ratio of maximum 
radial surface velocity to the maximum particle speed. A 
smaller ratio denotes a more adiabatic shape changes in 
relation to the particle speed. In analogy we define such 
index as a ratio 



>0?) 



Vf 



(15) 



of the root mean square radius speed in forward direction 
(opening angle i?) and the Fermi velocity. With opening 
angle 0.4 rad we obtain a maximum r\ = 0.12 for figure 
|l| and rj = 0.30 for fi gure 0. This shows that we are 
essentially still in the adiabatic regime described in ||. 

The nonlinear behavior described so far already doc- 
uments that we are in a regime of large amplitude os- 
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